Origin of the singular Bethe ansatz solutions for the Heisenberg XX Z spin chain 
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We investigate symmetry properties of the Bethe ansatz wave functions for the Heisenberg XX Z 
spin chain. The XXZ Hamiltonian commutes simultaneously with the shift operator T and the 
lattice inversion operator V in the space of = ±1 with the eigenvalue of T. We show that the 
Bethe ansatz solutions with normalizable wave functions cannot be the eigenstates of T and V with 
quantum number T) = (±1, =1=1) where T is the eigenvalue of V. Therefore the Bethe ansatz 
wave functions should be singular for nondegenerate eigenstates of the Hamiltonian with quantum 
number T) = (±l,=pl). It is also shown that such states exist in any nontrivial down-spin 
number sector and that the number of them diverges exponentially with the chain length. 
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I. INTRODUCTION 

The Bethe ansatz method has been applied to var- 
ious problems in condensed matter physics. The Heisen- 
berg XXZ spin chain, the Hubbard model, the Kondo 
model, and many other models are solved by the Bethe 
ansatz or its variants ||]. The Heisenberg XXZ spin 
chain describes various physical systems including inter- 
acting fermion systems. Also, it is related to the trans- 
fer matrix of the six- vertex model the asymmet- 
ric exclusion process, and a Kardar-Parisi-Zhang-type 
growth model For the XXZ Hamiltonian, the Bethe 
ansatz expresses energy eigenvalues and eigenstates as 
functions of quasi-particle momenta which should sat- 
isfy a set of algebraic equations, the Bethe ansatz equa- 
tions (BAE). Various numerical and analytic methods 
have been developed to solve the BAE, from which the 
ground state energy and low-lying excitation spectrum 
are calculated ||5|-|l0(| . Combined with a finite-size-scaling 
theory, they provide useful tools to study equilibrium and 
nonequilibrium critical phenomena. 

When applying the Bethe ansatz method to the 
Heisenberg XXZ spin chain, one may doubt whether the 
Bethe ansatz produces the complete set of eigenstates. 
The completeness of the Bethe ansatz is proved for an in- 
homogeneous generalization of the XXZ chain threaded 
by the Aharonov-Bohm flux Jll| , p^ . It implies that the 
Bethe ansatz for the XXZ Hamiltonian is also complete 
in a general sense. One may obtain the eigenstate of the 
XXZ Hamiltonian from that of the generalized Hamilto- 
nian by taking a limit where the inhomogeneity and the 
flux vanish. 

On the other hand, it has also been investigated for a 
long time whether the Bethe ansatz produces all eigen- 
states of the XXZ Hamiltonian without the help of the 
limiting procedure. There have been the arguments for 
the completeness assuming a string conjecture. For the 
XXX chain, which has isotropic spin-spin interaction, 
Takahashi showed that the number of the states con- 
structed from the Bethe ansatz using the string conjec- 
ture is equal to the total number of the spin states. The 



string conjecture was also adopted to prove the complete- 
ness of the Bethe ansatz for a generalized XXZ chain p^ ] 
and the Hubbard model Later, the string conjecture 
was shown to be invalid by Essler et al. |p6[|, who ana- 
lyzed the algebraic Bethe ansatz equations in the two 
down-spin sector. They found that some solutions ex- 
pected from the string conjecture are missing and that 
there exist solutions violating the string conjecture. Nev- 
ertheless, they concluded that the Bethe ansatz is com- 
plete at least in the two down-spin sector because the 
numbers of the former and the latter are the same. 

However one of the solutions found in Ref. [ p^ is singu- 
lar, the corresponding Bethe ansatz wave function is not 
well-defined. In fact the state |*) = + 1)> 

which is an eigenstate of the XXZ Hamiltonian, is not 
produced from the regular (non-singular) solutions of the 
Bethe ansatz (see below for the definition of the state 
1^, / -|- 1)). Quasi-particle momenta corresponding to the 
state are divergent |]l|,^18| and the Bethe ansatz wave 
function is ill-defined. An eigenstate that is not produced 
from the regular solutions of the Bethe ansatz will be re- 
ferred to as a singular state. Siddharthan reported 
that there also exist the singular states in the three down- 
spin sector. It was claimed that a symmetry property 
might be important since the singular states have a def- 
inite symmetry property. 

The singular state does not imply the incompleteness 
of the Bethe ansatz for the XXZ Hamiltonian. One may 
obtain the wave functions for the singular states of the 
XXZ Hamiltonian from those of the generalized Hamil- 
tonian by taking the limit where the inhomogeneity and 
the flux vanish (see, for example, Ref. ||l^]). However it 
should be examined whether the singular states are rele- 
vant to the involved physical quantities of the Heisenberg 
XXZ spin chain since, if so, one must first consider the 
generalized model and take the suitable limit . 
Therefore, it is helpful to answer the question why there 
appear the singular states and how many there are in 
general down-spin number sectors. In this paper, we in- 
vestigate the symmetry properties of the XXZ Hamilto- 
nian and the Bethe ansatz wave functions to show that 
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the regular solution of the Bethe ansatz with well-defined 
wave functions does not contain a certain class of eigen- 
states. Those states are the singular states. We also es- 
timate the number of the singular states approximately 
at each sector with fixed number of down spins and find 
that the number of them diverges exponentially with the 
chain length. 



II. BETHE STATES AND SYMMETRY 
OPERATIONS 

The Heisenberg XX Z spin chain of length N has the 
Hamiltonian 
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H 



4i: 



1=1 



(1) 



where af'^'^ are the Pauli matrices at site — 
{af ± iaf)/2. The periodic boundary condition is im- 
posed, crJv+i ~ cFi'^'^- We assume N is even unless stated 
otherwise. It reduces to the ferromagnetic (antiferromag- 
netic) XXX chain in the case of A = 1 (A = —1). 
Since the Hamiltonian commutes with the magnetization 
M = , one can work in the Q down-spin sector in 

the diagonal basis of cr^. A state in the Q down-spin sec- 
tor is spanned by nCq state vectors \xi, X2, ■ ■ ■ , xq) with 
1 < Xi < N {xi < Xj for i < j) denoting the position of 
the ith down spin. 

The eigenstate of H can be written as j^*) — 
T,xi< - <XQ '/'(^^i, • ■ • ,2:Q)|a:;i, . . .,xq). According to the 
Bethe ansatz, the wave function is given in terms of quasi- 
particle momenta km = 1, ■ • ■ , Q) by 



1p{xi,...,XQ) 



A{P) exp j i ^ kp^Xr, 

P \ m=l 



(2) 



where the sum is over all permutations P of integers 
{1,2,..., Q}. The quasi-particle momenta are complex 
numbers with — tt < Re k < n. The amplitudes ^'s for 
different permutations P and P' , which are identical ex- 
cept for two neighboring integers such that Pj = -Pj+i = 
TO, Pj+i ^ Pj — I, and Pi — P- for i ^ + are related 
as 



MP') 



gi{k^+k,) _ 2Ae'''' + 1 



(3) 



The periodic boundary condition implies the Bethe 
ansatz equations (BAE) for the quasi-particle momenta: 
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j.{k^+ki) _ 2Ae*'^'" -I- 1 



(4) 



The state constructed from the Bethe ansatz is called the 
Bethe state. It has the energy E — ~ cosfc^)- 



In the special case of A = 0, where the BAE becomes 
gik^N _ (|_2)Q-i^ Bethe ansatz produces all eigen- 
states of H. A solution is obtained by selecting Q differ- 
ent values among the N values of (7^)x integers (half- 
integers) for odd (even) Q. There are nCq different ways 
for each Q, which is equal to the possible spin states. 
Naively one may think that the solutions at A = evolve 
continuously as A is turned on. However it will turn out 
to be false. 

The XXZ Hamiltonian under the periodic boundary 
condition commutes with the shift operator T, which 
shifts the position of each down spin to the left by one 
unit. By construction, the Bethe state is the eigenstate 
of T with the eigenvalue ft = [11™=! e'M- Since is 
the identity operator, fl^ should be unity. The Hamil- 
tonian also commutes with the lattice inversion operator 

V defined by T^|{a;J) = \{yi}) with yi = N~xq-i+i + 1. 
Since is the identity operator, the eigenvalue of V, 
which will be denoted by T, takes the value of ±1. The 
inversion operator does not commute with the shift oper- 
ator in general. But it is easy to show that they commute 
with each other in the subspace oi H. — ±1; [T,V]\'^) = 
if r|vl') — ±1^). Then the eigenstate of H can be made as 
the simultaneous eigenstate of T and V in the subspace. 

The Bethe state with = ±1 is not the eigenstate of 

V necessarily. If a Bethe state with = ±1 with 
well-defined wave function is a simultaneous eigenstate 
of (V^l^f) = T|$)), then one can show that T should 
be equal to fl : 



T = n . 



(5) 



It is proved as follows. Suppose that is a Bethe state 
in the Q down-spin sector with the quasi-particle mo- 
menta {fci, ^2, • . . , fcg} satisfying Eq. (|j). The wave func- 
tion is given as in Eq. (|^). When one applies the lattice 
inversion operator V to the state |^), the wave function 

of the state = J2 i^i^i^ • ■ • j ^q)\^1j ■ • ■ j ^q) is given 

by 



iPixi, ...,xq)^ iJjiN -XQ + l,---,N-xi + l) 



= n'^AiP) exp 



i^(-fcp„)a 



, (6) 



where A{P) = A{PP) with a parity permutation P 
which maps {1, 2, • • • , Q} to {Q, •••,2,1}, and we used 
that [11™=! e*'^"] = Q and fl^ ~ 1. One can easily ver- 
ify that {—k,n} satisfies the BAE in Eq. (^) and that the 
ampUtude A{P) satisfies Eq. (^) with {km} replaced by 
{—km}. This shows that V\'$) is a Bethe state with the 
quasi-particle momenta {— fci, — ^2, . . . , — kg}. Therefore 
a Bethe state could be an eigenstate of V only when 
the two sets {/ci, fc2, . . . , fcg} and {— fci, — ^2, . . . , — fcg} 
are identical. 

We rewrite the wave function in Eq. (|^), introducing 
the permutation P defined by the relation —ki = k~ , as 
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^ ^ A{PPP) exp 



(7) 



Comparing it with Eq. (|^), one can find T from the ra- 
tio of A{P) and flA{PPP) for a certain permutation P, 
e.g., the identity permutation /. If the amphtude A is 
nonzero and finite, T is given by 



T = nA{PP)/A{I). 



(8) 



Note that the set {km} is equal to {—km}- On one 
hand, unless both and tt are present in the set {km}, 
one can rearrange the momenta in such a way that the 
permutation P and the parity permutation P are the 
same. Then one has A{PP) = A{P'^) = A{I) and hence 
T = il. On the other hand, if both and tt are present 
in {km}, Q should be even and one can rearrange the 
momenta in the following way: 

{ki, .. ., fcQ/2-l, ^Q/2+l, -kQ/2-i, ■ • ■ , -ki} , 

with kQ/2 = and fcQ/2+1 = With this arrangement, 
the permutation PP maps Q/2 to Q /2 + 1, (5/2 + 1 to 
g/2, and j ^ Q/2, Q/2 + 1 to j. Then, using Eq. (|), 
the ratio between A{I) and A{PP) is given by 



A{PP) 



= 1. (9) 



Therefore we find that = Jlj^') if |*) is an eigen- 

state of V. This proves Eq. 



III. PERTURBATIVE CONSTRUCTION OF 
SINGULAR STATES 

Equation (||) of previous section shows that the Bethe 
state with well-defined wave function cannot be the 
simultaneous eigenstate of T and V with (O, T) ~ 
(±1,=f1). But there exist many eigenstates of H with 
such symmetry. Suppose an eigenstate of H possesses the 
quantum numbers (ri,T) = (1,-1). The state may be 
degenerate (having the same energy) with another state 
which has T = -1-1 in the same sector. In the degener- 
ate case, the Bethe state may be a linear superposition 
of the degenerate pair, which does not lead to the contra- 
diction to Eq. (^. On the other hand, if the eigenstate 
does not have the degenerate partner with T = +1 in 
the same ft sector, such a state cannot be a Bethe state, 
hence is absent in the regular Bethe ansatz solutions. In 
the remaining part, we show that there exist nondegen- 
erate states with (fi, T) = (±1, =f1) for nonzero A in the 
general Q down-spin sector. 

It is convenient to work with the second quan- 
tized form of the Hamiltonian. Under the Jordan- 
Wigner and Fourier transformations, defined as = 



a; exp[i7r^^-^\ ajoj] and ai = X^p^p^'^VV^, respec- 
tively, the XX Z Hamiltonian in the Q down-spin sector 
reads as H = Hq + Hi where Hq = — X]p(cosp — A) Up 
and 



P1+P2=P3+P4 



Here Up and aj, are anticommuting fermion operators 
with the bare momentum — vr < p < tt and Up = a^^ap is 
the number operator with "^pUp = Q. The momentum 
p takes the real value of {2tt/N)x integer (half- integer) 
for odd (even) Q. 

Consider a symmetric set Sq-2 which consists of 
{Q — 2) bare momenta satisfying J2peSQ-2 P — ^ (mod 
2n). A set S is symmetric if it contains —p for all p G S. 
For a given Sq-2, R is defined as the set of all momenta 
in the interval — 7r/2 < p < 7r/2 except p' or n — p' (mod 
27r) for p' G S'q-2- Note that R is also a symmetric set. 
Then, at A = 0, the following states 



\pa; S, 



Q-2, 



= t t 
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AM 



(10) 



with Pa & R are degenerate eigenstates of H in the Q 
down-spin sector having Vl = —1. Here |0) denotes the 
vacuum, i.e., ap|0) = for all p. They are not the eigen- 
states of V . One can easily show that 



V\pa\SQ-2 



-VL\ -pa, Sq-2 



(11) 



Figure [T] illustrates an example of the set Sq-2 and cor- 
responding R with the states in Eq. ( p^ for = 20 and 
Q= 10. 

As A is turned on. Hi generates the overlap between 
the degenerate states, which may cause the degeneracy 
to split. To see the degeneracy splitting, we treat Hi 
as a perturbation and apply the degenerate perturba- 
tion theory. In this scheme, the leading order correction 
to the energy is given by the eigenvalue of the pertur- 
bation matrix (iJi)p^^p^ = {pa\ SQ-2\Hi\pp; Sq-2) with 
PcnPii G R- The eigenstates are given by the eigenvectors 
of the perturbation matrix. A straightforward calcula- 
tion shows that the matrix element is given by 

4A 

(-ffl)p„,p^ ^ - — {cOSPa){cOSpi3) , 

where we neglect a constant diagonal element. It has the 
nondegenerate eigenvector 



\^\Sq-2)= X! C0Spa\pa\SQ-2) 



(12) 



with the eigenvalue —(4 A/TV) ^nCOS^Pa- The other 
states remain to be degenerate with their common eigen- 
value 0. If one applies V to the state Sq-2) and uses 
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Eq. (p, l|), one can see that the nondegenerate state is the 
eigenstate of V with T = +1 while ^l = —1. Therefore, 
it cannot be produced from the regular solution of the 
Bethe ansatz by Eq. (||). 

The remaining degenerate states may split further in 
the higher order perturbations. So the number of the 
nondegenerate states with {ft,T) = (— 1,+1) is equal to 
or greater than the number of different realizations of 
the set Sq-2- For even Q, 5q_2 consists of {Q — 2)/2 
values of p selected randomly in the interval < p < tt 
and their negative partners. For odd Q, Sq-2 consists 
of p = 0, (Q — 3)/2 p's randomly-selected in the inter- 
val < p < TT, and their negative partners. So, there 
are at least n/2C(q^2)/2 (7v/2C'(q-3)/2) states which have 
(ri, T) = (—1,-1-1) and they are those states which can- 
not be produced from the Bethe ansatz in the even (odd) 
Q down-spin sector. If we sum over all Q, the number of 
such states is of the order of 2^/^. 

In the two down-spin sector, the perturbation calcu- 
lation becomes exact since the states in Eq. ( |l0| ) with 
the empty set S'q-2 span the whole two down-spin sector 
with $7 = — 1. So the state 

Yl «p4-p|0) 
-f <p<f 

is the eigenstate of H at all A, which cannot be pro- 
duced from the Bethe ansatz. If one performs the in- 
verse Jordan- Wigner and Fourier transformations, the 
state can be written as ^[{—^Ylhi + !)■ This is the 
example of the singular states mentioned in |l7| . 
Consider also the states 



IV. NUMERICAL RESULTS 

The state counting discussed in the previous section is 
not exact. States from different S sets may be degenerate 
at A = 0. In that case perturbation calculations should 
be carried out in an enlarged space and it is not guar- 
anteed that each S set generates at least one singular 
state. And the singular states with (17, T) = (±1,=f1) 
may become degenerate accidentally with a state with 
(f2, T) — (±l,±l) for a particular value of A. Then the 
Bethe ansatz may produce the linear superpositions of 
them. We believe that those effects do not affect the 
leading order of magnitude of the total number of the 
singular states. 

We checked this by diagonalizing the XX Z Hamilto- 
nian exactly with A'^ = 6,8,...,16. The symmetry prop- 
erty of each eigenstate is checked and the nondegenerate 
states with (J7, T) = (+1, —I) and (—1, +1) are identified 
as the singular states. Figure || shows the energy level, 
as a function of A, of each state with O = ±1 for N = 8 
in the four down-spin sector. The singular states are 
denoted by dotted lines. Note that each singular state 
evolves from the degenerate states at A = as the de- 
generacy is lifted for nonzero A . At A = ± i , many sin- 
gular states become degenerate with non-singular states 
and they are not counted as the singular states. In Ta- 
ble 1, we show the total number of the singular states at 
A = 0.2,0.5, and 1.0. Due to the accidental degeneracy 
the number varies with A. It is confirmed that the num- 
ber increases approximately as 2^^^ at the three values 
of A. 
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with Pa,Pi3 G R, constructed from a symmetric set 5q_4 
consisting of (Q — 4) momenta satisfying J2peSQ^iP ~ 
(mod 2tt) and the set R of all momenta in the interval 
— I < p < § except p' or n — p' (mod 2tt) for p' G Sq-4. 
They are degenerate eigenstates of at A = with 
n = +1. The same perturbation analysis will show that 
the degeneracy splits for finite A and there appear non- 
degenerate states with T = — 1. Those states cannot be 
produced by the Bethe ansatz from Eq. (^) either. There 
are n/2C{q-4)/2 {n/2C{q^5)/2) different ways in choosing 
the set Sq-4 for even (odd) Q. For each set, there exists 
at least one state with (il, T) = (-1-1,-1). So the total 
number of such states is again of order of 2^/^. 

The same analysis can be applied to the states from the 
symmetric sets Sq-q, Sqs, . . ., with X^peS-P — (mod 
27r) and symmetric sets S'g-i, S'q_3, . . . containing p = tt 
and satisfying J2pes ~ 2tt). For each case, there 

are at least a number of order of 2^/^ singular states. 



The singular behaviors of some Bethe ansatz solutions 
have been noticed for the Heisenberg XX Z model (see, 
e.g., p],p|jT7[ ) . So the generalizations to the inhomoge- 
neous XX Z model with Aharonov-Bohm flux were con- 
sidered and their Bethe ansatz solutions were shown to be 
complete [|ll|,0 . In those cases the Hamiltonian does not 
possess the lattice inversion symmetry and the present 
analysis is not applicable. 

The singular states include low-lying states in the half- 
fiUing case, Q = N/2. Figures 1 (b) and (c) illustrate 
degenerate first excited states at A = in the = -1-1 
sector. They are not the eigenstates of the inversion op- 
erator V. For nonzero A the degeneracy is removed, and 
the states evolve into the eigenstates of V with T = ±1. 
Therefore the Bethe ansatz fails to produce one of them 
with T = — 1. Those states are important in identifying 
the operator content of the XX Z model in the context of 
the conformal field theory (see j^] and references therein). 
They have the scaling dimension Xn.m with n = and 
m = ±1. To recover the whole operator content from the 
Bethe ansatz, one should introduce the Aharonov-Bohm 
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flux and take the limit where the flux vanishes as done 



So far we have assumed N is even. For odd N, one 
cannot find a pair of bare momenta p and p' satisfying 
p + p' = TT and the previous analysis cannot be applied. 
So the above arguments are not valid for odd N, though 
Eq. ^ holds good. We have performed the same numer- 
ical analysis for odd TV < 15 and found that the nonde- 
generate states with (fi, T) = (±1, =f1) do not exist. 

In summary, we have shown that the symmetry prop- 
erties of the Hamiltonian and the Bethe wave functions 
prevent the regular solutions of the Bethe ansatz from 
containing some eigenstates. For the fixed number of 
down spins, the XXZ Hamiltonian can be simultane- 
ously diagonalized with the shift operator T and the 
lattice inversion operator V in the subspace where the 
eigenvalue of T, Q, is ±1, and each simultaneous eigen- 
state has one of the following four pairs of eigenvalues, 
{n,T) = (1,1), (-1,-1), (1,-1), (-1,1). But there is 
no regular Bethe ansatz wave function with (fi, T) = 
(1, —1), (—1, 1), which leads to the appearance of the sin- 
gular Bethe states. We estimated the number of the sin- 
gular states as of 0(2^/^) for even TV. The number of 
the singular states is exponentially large though the frac- 
tion to the total number of the states, 2^, is vanishingly 
small. Moreover, the singular states include low-lying 
states in the half-filling case where the number of down 
spins is N/2. 
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TABLE L The number of the nondegenerate states with 
(n,T) = (±l,Tl). 
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A = 0.2 



A = 0.5 



A = 1.0 



8 

10 
12 
14 
16 



12 
25 
64 
152 

339 



7 

15 
29 
61 

131 



12 
20 
59 
152 

319 
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(a) (b) 
1. An illustration of sets 



(c) 

FIG. 1. An illustration of sets Sq-2 and R, and a state 
in Eq. with = 20 and Q = 10. The bare momenta 
p are denoted by equi-spaced circles in a complex z = e'*' 
plane. The horizontal (vertical) dashed lines stand for the 
real (imaginary) z axis, (a) The elements of the set Sq-2 (R) 
are denoted by black (gray) circles. In (b) and (c) black circles 
denote the occupied momenta in Eq. (id). 
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FIG. 2. Energy vs. A with = ±1 of the XXZ Hamilto- 
nian in the four down-spin sector with the chain length N = 8. 
The dotted lines represent the singular states, nondegenerate 
eigenstates with quantum number (f2, T) — (±1, =Fl)- 
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